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TeV-PeV Neutrinos from Giant Flares of Magnetars and the Case 

of SGR 1806-20 

Kunihito loka} Soebur Razzaque^ Shiho Kobayashi^’^, and Peter Meszaros^’^ 

ABSTRACT 

We estimate the high energy nentrino flux from the giant flare of SGR 1806- 
20 on December 27, 2004, which irradiated Earth with a gamma-ray flux ~ 10^ 
times larger than the most luminous gamma-ray bursts (GRBs) ever detected. 

The Antarctic Gherenkov neutrino detector AMANDA was on-line during the 
flare, and may either have detected high energy neutrinos for the first time from 
a cosmic point source, or put constraints on the flare mechanism of magnetars. 

If TeV neutrinos are detected, one would expect also detectable EeV cosmic rays 
and possibly TeV gamma-ray emission in coincidence. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory — 
stars: individual (SGR1806-20) — stars: neutron 


1. Introduction 

The giant flare of SGR 1806-20 on December 27, 2004 was the brightest cosmic transient 
to date, emitting gamma-rays for ~ 0.1 sec with flux ~ 10 erg s“^ cm“^ (Terasawa et al. 
2005; Hurley et al. 2005; Mazets et al. 2005; Palmer et al. 2005). If a comparable energy 
was emitted as high energy neutrinos, they should have been detected for the first time 
by current neutrino observatories such as AMANDA (Ahrens et al. 2002). In this Letter 
we calculate the expected high energy neutrino flux from giant flares in Soft Gamma-ray 
Repeaters (SGRs) such as SGR 1806-20, and argue that the neutrino flux is indeed either 
detectable, or its absence provides important constraints on the flare mechanism (see Zhang 
et al. (2003) for high energy neutrinos from quiescent magnetars). 
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SGRs are a type of extreme X-ray pulsars, repeatedly emitting ~ 0.1 sec bursts of soft 
gamma-rays. Giant flares are more energetic events, which have been recorded from three 
of four known SGRs with the December 27 event being the third one. SGRs are most likely 
magnetars, i.e., highly magnetized neutron stars (Thompson & Duncan 1995, 2001). In 
this model giant flares result from a global magnetic rearrangement of the crust or even the 
entire interior (e.g., loka 2001). 

Giant flares have many similarities to cosmological gamma-ray bursts (GRBs). Long 
duration GRBs are thought to arise from relativistic jets interacting with themselves for 
prompt GRBs, and subsequently with a circumburst medium for the longer wavelength 
afterglows (e.g., Zhang & Meszaros 2004). In a similar manner, the detected radio afterglows 
of SGRs imply the presence of relativistic outflows (Frail, Kulkarni & Bloom 1999; Gameron 
et ah 2005; Gaensler et ah 2005), and the huge flare luminosities also lead to relativistic 
fireballs (Huang, Dai & Lu 1998; Thompson & Duncan 2001; Nakar, Piran & Sari 2005). 
In particular, the minimum energy of the radio afterglow is larger than the kinetic energy of 
pairs that survive annihilation, implying the presence of baryons in hreballs like GRBs 
(Nakar, Piran & Sari 2005). Relativistic baryons are likely to cause shocks, leading to Fermi 
accelerated protons, and to high energy neutrinos via py (Waxman & Bahcall 1997) and 
pp interactions (Paczyhski & Xu 1994). The accelerated protons in GRBs may also explain 
the observed ultrahigh energy cosmic rays (Waxman 1995; Vietri 1995; loka, Kobayashi, & 
Meszaros 2004). 

In contrast to GRBs, however, the baryon load in giant SGR flares is less constrained (see 
§ 2), mainly because the flare spectrum may be thermal (Hurley et al. 2005) or nonthermal 
(Mazets et al. 2005; Palmer et al. 2005). Since the neutrino fluxes depend on the baryon 
load, neutrinos can be a probe of the baryon load as well as the flare mechanism. In § 2 we 
discuss typical hreball models for SGR 1806-20 and in § 3 we estimate the expected high 
energy neutrino fluxes for these models. The implications are discussed in § 4. 


2. Typical fireball models for SGR 1806-20 

The giant flare of SGR 1806-20 on December 27, 2004 radiated a gamma-ray energy 
F'.y ~ 3 X 10"^® 77^^46,5 erg during a time of to ~ 0.1to,-i sec (Terasawa et al. 2005).^ The 
total luminosity was Lq ^ 7/^7 ~ 3 X 10"^'^Lo, 47.5 erg s ^ for a conversion efficiency of 
total energy into gamma-ray. If such energy is released near a neutron star (radius ro ~ 


^We adopt d = lOdi kpc for the distance to SGR 1806-20 though it is controversial (Cameron et al. 2005; 
Corbel & Eikenberry 2004). 
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10 ®ro,6 cm), pair production creates an optically thick fireball with an initial temperature 
(Paczyhski 1986; Goodman 1986) 
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where a = 7r^/c^/15h^c^ = 7.6 x 10“^^ erg cm“^ K“^ is the radiation density constant. As 
the hreball expands under its own pressure, the Lorentz factor increases as F oc r and the 
comoving temperature drops as T oc r~^. The subsequent evolution depends on the baryon 
load parametrized by a dimensionless entropy r] = Lq/Mc^ (Shemi & Piran 1990). If the 
hreball is baryon-rich, r] < r]^, where f]^ is a critical entropy (Meszaros & Rees 2000) 
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a photosphere appears in the coasting phase and almost all the energy goes into the kinetic 
luminosity of the outhow Akin ~ Lq, while a photosphere appears in the acceleration phase 
if ?7 > ? 7 *. We can derive 


Tf = min[? 7 ,? 7 *], rph/rg = max[r/*(7//?7*) ^, r]^{r]/r]^) 

Lkin/Lo = min[l,? 7 */? 7 ], Lph/To = min[(? 7 /? 7 *)®/^, 1], Tph/Tg = min[(? 7 /? 7 *)®/^ 1], (3) 


where the hrst (second) value in the bracket is for r] < [f] > i/*), P/ is the hnal Lorentz 

factor, a thermal photosphere at radius rph has a luminosity Lph with an observed tempera¬ 
ture Tph, and we neglect hnite shell effects for simplicity (Meszaros et al. 2002). The above 
relations hold provided the hreball is not too baryon-poor, ?/ > ~ lO^L^y^grg g^ (Meszaros 

& Rees 2000). 

Internal shocks in a variable outhow with P = 10^r2 are expected to occur at radii 


Vg — 2P^cAt ~ 6 X lO^^P^At.i cm (4) 

where At = 10“^At_i s is the variability timescale. The minimum value is At r'U ro/vA rsj 
3 X 10“^ s since the Alfven velocity in the magnetosphere is va ~ c, which is consistent with 
the initial rise time < 0.3 ms (Terasawa et al. 2005; Palmer et al. 2005). If shocks occur 
above the photosphere, i.e., rj > rjs = IOLq nonthermal radiation is produced. 
The nonthermal shock luminosity Lg = ■Cs-^kin dominates the photospheric luminosity Lph if 
7] < 1007/*^2Cs^*) where is the conversion efficiency of kinetic energy into photons. 

Two typical scenarios are possible for the December 27 hare. The hrst one is a baryon- 
poor scenario, e.g. r] ~ 10"^ > ?/*, where the photosphere is in the acceleration phase. Most 
of the energy is radiated as photospheric emission, which may explain the thermal spectrum 
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with temperature Ty ~ 175 ± 25 keV observed by Hurley et al. (2005). The remaining 
kinetic energy is ~ 10 ^ times smaller than the radiation energy, which is also implied 
by the radio afterglow of SGR 1806-20, if we use typical parameters inferred from GRB 
afterglow httings 0.15 0.01 and n ~ 1 cm ^ (Nakar, Piran & Sari 2005; Wang et 

al. 2005). 

The second is a baryon-rich scenario, e.g. with t] ^ 10 < r]^,. The observed thermal 
radiation can be explained by a thermalized shock luminosity Lg = ^sL]nni if shocks 
occurred right below the photosphere (Rees & Meszaros 2004). The remaining kinetic energy 
is (1 - ^.)/6 10.^3 times larger than the radiation energy, so that the observed radio 

afterglow may reqnire a jet conhguration {9 ;< 0 . 1 ) or atypical model parameters (high 
ambient density n and low ,^e or ^b)- However a jet may be implied by the elliptical image 
and polarization of the radio afterglow (Gaensler et al. 2005) as well as the light curve 
of the giant flare that is well htted by emission from a relativistic jet (Yamazaki et al. 
2005). Atypical parameters may be also snggested by the rapid decay of the radio afterglow 
(Gameron et al. 2005). Thus this model may also be viable. 

The actual spectrum of the November 27 flare, however, may be nonthermal (Palmer et 
al. 2005; Mazets et al. 2005). At least a portion of the giant flare may be nonthermal, since 
Hurley et al. (2005) determined the spectrum with low time resolution. A previous giant 
flare in SGR 0526-66 may also have had a nonthermal spectrum (Fenimore, Klebesadel, 
& Laros 1996). If giant flares are observed as short GRBs, their spectra are also likely 
nonthermal (Fenimore, Klebesadel, & Laros 1996; Nakar et al. 2005). A nonthermal flare 
component may arise in a baryon-rich model rj > 10 , since internal shocks extending above 
the photosphere naturally produce nonthermal emission. The characteristic synchrotron 
freqnency = Tq^hgR/meC in internal shocks can be estimated as 

~ 000eji\TLlll,^,V-,^At-_\ keV, (5) 

where a fraction of the internal energy goes into electrons ( 7 ^ ~ ^eT^p/^e) and a fraction 
goes into the magnetic held, 47 rr^cF^R^/ 87 r = Note that this synchrotron 

emission is not thermalized above the photosphere. 


3. Proton interactions and nentrinos 

According to § 2, we choose the particular parameters for the two typical models as 
baryon — poor (BP) : 7 ~ 10 ^, Akin ~ 10~^L^ ~ 10 "^^'^ erg s“^, P ~ 100 , At ~ 10 “^ s, 
baryon — rich (BR) : 7 ~ 10 , Akin ~ lOA^ ~ 10 "^®'^ erg s“^, F ~ 10 , At ~ 10 “^ s. 
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We set At in the baryon-poor model to maximize neutrino events, which turn out to be 
undetectable. In the baryon-poor (baryon-rich) model internal shocks take place outside 
(inside) the photospheric radius (rxh ~ !)• Since the internal shocks are mildly relativistic, 
protons are expected to be accelerated to a power law distribution, dup/dep oc (Waxman 
& Bahcall 1997; Waxman 1995; Vietri 1995). High energy protons may then interact 
with photons (py) or with other cold protons (pp) to produce neutrinos mostly through 
pion decays. The maximum shocked proton energy in the comoving equipartition magnetic 
held: B' = -^^’nn'p'mpC^^Bii) where n'^ is the comoving baryon number density, is limited by 
the system size and synchrotron losses. By equating the acceleration time ~ e'p/cqB' 
to the shorter of the comoving time ~ Ts/Tc and the synchrotron cooling time ~ 
GTirripC^/aTmlepB"^ we hnd the maximum proton energy in the lab-frame as 


f- r\j 

^p,max 


4 X 1016 (eB,- 26 ,-iW 5 . 5 )'/'r 2 -i eV (BP) 

7 X IQis (eB,_26,-ii^kin,48.5)-'/"r^/'AfL/i' eV (BR), 


( 6 ) 


where a fraction ~ 0.1.^i^_i of the kinetic energy goes into the internal energy. Note that 
is smaller (greater) than in the baryon-poor (baryon rich) model. 

The comoving baryon density in the two models is 


= 


^kin 


2 X 10iiLkin,45.5r2'^At_4 cm ^ (BP) 
I’ “ 47rr2p2mpc3 ' \ 2 x lOi^Lkinys.srr^Atl^ cm'^ (BR). 


(7) 


The accelerated proton hux that would be measured at Earth, if they were to reach us in a 
straight line before converting to neutrinos is 




20 (ep/GeV)-2e,,_iLkin,45.5C^r' GeV"! cm’^ s'l 


(BP) 


2 X 10^ (ep/GeV)-2^i,_iLkin,48.5t^r^ GeV"! cm-2 s"! (BR). 


( 8 ) 


A fraction of this proton flux will be converted to neutrinos depending on the opacity of py 
and pp interactions in the hreball. 

The photospheric thermal radiation bathes the ejecta at the internal shock region (which 
also radiates), so that accelerated protons interact with photons (using the observed T^) of 
comoving energy 


4 ~ r,/r ~ P ry kev (bp) 

’ \ 20 Pf' keV (BR) 

and can produce pions if the observed proton energy is 


(9) 
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( 10 ) 
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which is below the maximum available proton energy in both the models. The density of 
these thermal photons at the shocks is^ 


/ / 7 X 10 ^* 1 .^, 47 , 56 .^ 5 312^ At _4 cm ^ 

^ dTrr^r^ce; ~ 1 7 x L^,^r.5e:;i3rfAtzl cm"* 

The corresponding optical depth to py interactions is then 

(Tp^n'^r, f 2L^,47.5S,k3r7^7it:^ (BP) 
~ T ~ \ 20L.,,47.5e-k3rr"At:} (BR), 


(BP) 

(BR). 


( 11 ) 


( 12 ) 


where apy ~ 5 x 10 cm^ is the cross-section at the A resonance. 


Protons lose ~ 20% of their energy at each pj interaction, dominated by the A reso¬ 
nance. Approximately half of the pions are charged and decay into high energy neutrinos 
77/i+ -|- ^ e"*" -|- z/g + -|- z/^, with the energy distributed roughly equally among the 

decay products. Thus the neutrino energy is ~ 5% of the proton energy. From equation 
( 10 ), we find the typical neutrino energy expected from the pj interactions with thermal 
photons as 

O.OP^GeV^ f8xl05p2e;]3GeV (BP) 

206., ^ \ 8 X 10 *p 26 -^ 5.3 GeV (BR). ^ ^ 

The corresponding monoenergetic neutrino flux at Earth (equal for and z/g after oscil¬ 

lations in vacuum where Vp created from tt’*' decay is transformed to Vr) can be found from 
equation ( 8 ) as 


I . .0.2 

r 7 X 10 - 1 * 6,-iPkin,45.5e^,5.3^7''r2”" GeV' cm-2 s-1 (BP) 
\ 7 X 10-6 6,-iPkin,48.5e^,5.3^^r'r7" GeV' cm-2 s-1 (BR), 


(14) 


where in both the cases min(l,rp.y) = 1 from equation (12). Note that the synchrotron and 
inverse Gompton losses of tt’*' and pA are negligible for our models. 


In addition to py interactions, shock accelerated protons may also undergo pp interac¬ 
tions with cold protons in the ejecta and produce tt^. The opacity for pp interactions is 
about Tpp ~ appUpfs/T ~ 0.6 Lkinys.sPr^l^^^i ~ O.lrTh for an average pp cross-section of 
<^pp ~ 6 X 10-^6 in the TeV-PeV energy range. We calculate the neutrino flux from 
decays as (Razzaque, Meszaros, & Waxman 2003) 


= min(l,rpp) 



•Fp (6p) dcp, 


( 15 ) 


"‘Below the photosphere the photon density could be larger than this estimate. 
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where the neutrino multiplicity from pp interactions in units of GeV ^ is given by 
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(16) 


7cm is the Lorentz factor of the pp center of mass in the lab-frame and 0 is a step function. 
Note that the flux is given for muon (anti muon) neutrinos. Electron and tan (and their 
anti) neutrino fluxes would be the same on Earth. We may £t the pp muon neutrino flux for 


'pp 


< 1 by 




330 + 30 hi (^)] (^)'" «..-iiL.is. 5 rfr"rr'Arl GeV-'cm-=s-‘, (17) 


For the baryon-poor model with optically thin internal shocks, an additional neutrino 
component may arise due to p'y interactions with non-thermal synchrotron photons, but this 
latter flux component is undetectably low and is ignored. Next we calculate the expected 
neutrino events at Earth. 


4. Neutrino events in AMANDA 

SGR1806-20 is located in the southern sky at a declination of —20°. We have calculated 
the probability to detect muon neutrinos with the AMANDA detector, located at the South 
pole at a depth of 1 km, by using a code which propagates neutrinos through Earth and 
calculates the interaction rate in ice near the vicinity of the detector (Razzaque, Meszaros, 
& Waxman 2004). The resulting probability may be fitted in the TeV-EeV energy range 
with a broken power-law as 

P(e^) = 7 X 10-^(e^/10^-®GeV)^, (18) 

where f3 = 1.35 for < 10"^-® GeV while f3 = 0.55 for > 10"^'^ GeV. Using a geometrical 
detector area of Adet = 0.03 km^ and the flare duration = 0.1 s, the number of muon 
events from py neutrinos (z/^) of energy given in equation (13) is 

^ r 0.007 min(l,rp^)^i _iLkin,45.5C?r^(e^,Lri)^"^Adet-1.5^0-1 (BP) 

\ 20 min(l,rp..y)^i_iLkin,48.5'^r^(S,k3r?)^"^Adet-1.5^0-1 (BR). 

On the other hand, the number of muon events from pp neutrinos would be 

^PeV 

^^i,pp 2Adet^0 / Pi^^ii)^v,ppd(^u 

^TeV 

13 X 10“® .^j^_iLkin, 45.5'^! ^B2 ^Af_4Adet,-l.5^0, -1 (BP) 

\ 300 _lLkin,48.5'^l Bi Af_4Adet,-l.50,-1 (BR), 


rs_/ 


( 20 ) 



where we use equation (17) and note that Cherenkov detectors do not distinguish between 
and 9^ flavors. The difference of the arrival time between neutrinos and gamma-rays would 
be within ~ to ~ 0-1 s. Neutrinos can precede gamma-rays if internal shocks occur deeply 
inside the photosphere. 

In Figure 1, based on similar calculations, we show the parameter space where AMANDA 
and the future detector ICECUBE (Ajet = 1 km^) (Ahrens et ah 2004) can detect more 
than one muon event in the plane of the baryon load r] and the variability timescale At. 
(Note that background events are negligible.) Baryons are rich (poor) for r] < 10^ [r] > 10^) 
while the corresponding Lorentz factor is r ~ r/ (F 10^) in equation (3). We can see 
that TeV-PeV neutrinos may have been already detected by AMANDA if the giant flare 
is baryon-rich [r] < 30), while a nondetection would suggest a baryon-poor fireball. This 
offers the exciting prospect of gaining independent information about the baryon load or the 
bulk Lorentz factor of the fireball, the efficiency of proton injection and energy dissipation 
in shocks and the variability timescale associated with the flare trigger. Such parameters 
would constrain the energetics and the radiation mechanisms inferred from electromagnetic 
observations. 

If the neutrinos from this flare are detected by AMANDA, one would expect ICECUBE 
to be able to detect less energetic flares in this and other galactic SGRs. One can show 
that the baryon-rich model with a flare ~ 10“^ times smaller than that considered here can 
produce about one event in ICECUBE, and the rate of such flares is about ~ 1/10 yr. 

Since sedimentation due to gravity causes heavier elements to stratify down, the surface 
tends to consist of lighter elements (Alcock & Illarionov 1980). An absorption feature in 
bursts from SGR 1806-20 was interpreted as due to proton cyclotron lines (Ibrahim, Swank, 
& Parke 2003; Ho & Lai 2001). However the magnetar surface may not contain hydrogen 
since hydrogen could burn very fast (Chang, Arras, & Bildsten 2004). If a fireball contains 
some heavy nuclei, photo disintegration processes with thermal photons may substantially 
reduce the maximum energy of the nuclei (Puget, Stecker, & Bredekamp 1976). A detailed 
calculation requiring a Monte Carlo simulation is out of the scope of this Letter. Here we 
present qualitative remarks if a nucleus of arbitrary mass number A and charge Z survives 
photo disintegration while being accelerated in the internal shocks. Their maximum energy 
would be oc Z for and oc A?Z~^A ^ for in equation (6). The 

threshold energy of nuclei for the pion production is ~ A times equation (10). The pions 
decay into neutrinos and the typical neutrino energy is almost the same as equation (13). 
Since the neutrino energy is the same, the neutrino flux is the same as equation (14). (Note 
that the optical depth for ~ A pion production by nuclei is also the same as equation (12) 
because the photomeson cross-section is ~ ATp.y.) Therefore the number of muon events in 
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equation (19) does not change much. 

External shocks that produce the radio afterglows may also accompany neutrino emis¬ 
sions. A simple application of the GRB afterglow (Wang et al. 2005), however, shows that 
the typical synchrotron frequency is too low to make the py interactions. 

Neutrons produced by pj interactions can reach us in a straight line without decay if 
their energy is larger than ~ 10^® eV (e.g., loka, Kobayashi, & Meszaros 2004), and may be 
observed as coincident cosmic rays. Since ~ 10^® eV protons interact with ~ 40r^e“Jg eV 
flare photons, the py optical depth is about Tp^ ~ 8 x 10 “^L.y, 47 , 5 Atl}e “^8 if we assume a 
thermal spectrum. Then the number of cosmic ray events in Adet ~ 10^ km^ detectors such as 
AUGER (Abraham et al. 2004) may be ~ ep^pTp^Adetto ~ 

Neutral pions produced by py and pp interactions decay into two gamma-rays with a 
flux and energy comparable to neutrinos. These gamma-rays might be detected by Milagro 
(Atkins et al. 2003) if the flare had occured in the northern sky and if the gamma-rays 
escape the emission region without making pairs. 
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Fig. 1.— Parameter regions where AMANDA and ICECUBE can detect more than one 
muon event are shown in the plane of the baryon load tj and the variability timescale At. 
Baryons are rich (poor) for r] < 10^ [f] > 10^) while the corresponding Lorentz factor is 
10^) in equation (3). Solid (dotted) lines are for py (pp) neutrinos. Photospheric 
thermal emission dominates the nonthermal emission Lph > ■Cs-^kin on the right of the long 
dashed line, while internal shocks occur below the photosphere TTh > 1 (also leading to a 
thermal spectrum) on the left of the dashed line. Thus the flare spectrum is nonthermal 
in the shaded region. We used a normalization .^sTkin + Tph = = 3 x 10"^^ erg s“^ and 

adopted = 0.01, = 0.1 and = 0.1. 







